Mem.  S.A.It.  Vol.  81,404 
©SAIt  2010 


Memorie  deiia 


Jet  interactions  with  magnetized  clouds 

Preliminary  results  from  PIC  code  and  large-scale 
hydrodynamic  simulations  for  AGN  jets 

J.H.  Beall^’^’^,  J.  Guillory^,  and  D.V.  Rose^ 

'  Space  Sciences  Division,  Naval  Research  Laboratory,  Washington,  DC 
^  College  of  Sciences,  George  Mason  University,  Fairfax,  VA 
^  Voss  Scientific,  Albuquerque,  NM 
St.  Johns  College,  Annapolis,  MD,  e-mail:  beall@sjca.edu 


Abstract.  When  astrophysical  jets  interact  with  the  ambient  medium  through  which  they 
propagate,  they  lose  energy  and  entrain  (and  accelerate)  that  medium.  Recent  observations 
of  such  interactions,  including  time-lapsed  “movies”  of  both  AGN  and  microquasar  jets,  can 
be  used  to  inform  models  of  the  jet-ambient-medium  interactions.  In  this  paper,  we  discuss 
some  aspects  of  these  jets,  including  the  mechanisms  of  their  propagation,  their  constitution, 
and  the  non-linear  character  of  their  energy  loss  via  plasma  processes.  We  also  present  for 
the  first  time  PIC  code  simulations  to  show  momentum  transfer  via  caviton  formation  as 
a  result  of  plasma  processes.  An  illustration  of  the  microphysics  of  the  interaction  of  a 
proton-electron  beam  with  an  ambient  magnetic  field  is  also  presented. 
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1.  Introduction 

Accretion  onto  a  massive  black  hole  can  pro¬ 
vide  such  large  amounts  of  energy  to  power 
the  jets  in  astrophysical  sources.  Wu  et  al. 
(2002);  Bower  et  al.,  (1995);  and  Beall  et 
al.  (2003)  have  noted  that  the  energy  pro¬ 
duction  rate  based  on  accretion  is  dEjclt  ~ 
rfil I2){dml df)c^  where  rj  is  an  efficiency  factor 
usually  taken  to  be  ~  10%,  and  c  is  the  speed 
of  light.  We  note  that  the  efficiency  for  nuclear 
reactions  is  ~  1%,  which  makes  conventional 
stellar  processes  an  unlikely  energy  source  for 
AGN. 
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The  luminosity  available  from  accretion  at 
the  Schwarzschild  radius  is,  therefore,  L  - 
rjA.SxlQ^'^dmldt  in  ergs  where  dm/dt  is  in 
gms  or  L  —  ri3xl0‘^dMoldt  in  ergs  s“', 
where  dMoJdt  is  in  solar  masses  per  year.  If 
the  source  has  a  luminosity  of  10"^“*  ergs 
and  T]  ~  0.1,  we  infer  an  accretion  rate  of 
3xlQ-^Moyr-K 

If  the  jet  persists  for  at  least  a  time,  Tjet,  of 
order  the  light  travel  time  along  its  length,  then 
Tjet  ~  3x10^  years  for  a  100  kpc  jet.  This  yields 
a  total  kinetic  energy  release  of  at  least  9x10^^ 
ergs  over  the  ‘lifetime”  of  a  large  scale  jet.  A 
Mpc-scale  jet  traveling  at  0.1c  would  suggest 
an  energy  of  9x10^®  ergs. 
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2.  Modeling  the  jet  interaction  with 
the  ambient  medium 

Analyses  based  on  hydrodynamic  simulations 
have  demonstrated  a  number  of  interesting  ef¬ 
fects  originating  from  ram  pressure  and  the 
consequent,  turbulent  acceleration  of  the  ambi¬ 
ent  medium  (see,  e.g.,  Basson  and  Alexander, 
2002;  Zanni  et  al.  2005;  and  Krause  and 
Camenzind  2003). 

However,  these  hydrodynamic  approaches 
neglect  an  important  species  of  physics;  the 
microscopical  interactions  that  occur  because 
of  the  effects  of  particles  on  one  another  and 
of  particles  with  the  collective  effects  that  ac¬ 
company  a  fully  or  partially  ionized  ambient 
medium  (i.e.  a  plasma).  A  detailed  discussion 
of  these  effects  can  be  found,  variously,  in 
Scott  et  al.  1980,  Rose  et  al.  1984;  Rose  et  al. 
1987;  Beall  1990,  and  Beall  et  al.,  2003. 

For  this  analysis,  we  posit  a  relativistic 
jet  of  either  e*,  p  —  e^,  or  more  generally,  a 
charge-neutral,  hadron-e“  jet,  with  a  signifi¬ 
cantly  lower  density  than  the  ambient  medium. 
The  primary  energy  loss  mechanism  for  the 
electron-positron  jet  is  via  plasma  processes, 
as  Beall  (1990)  notes.  Kundt  (1987,  2001)  also 
discusses  the  propagation  of  electron-positron 
jets. 

While  the  physical  processes  in  the  plasma 
can  be  modeled  by  PIC  (Particle-in-Cell)  codes 
for  some  parameter  ranges,  astrophysical  ap¬ 
plications  of  the  PIC  code  are  not  possible  with 
current  or  foreseeable  computer  systems.  We 
therefore  model  these  plasma  processes  for  the 
astrophysical  regime  by  means  of  a  system  of 
coupled,  differential  equations  that  represent 
the  normalized  wave  energy  densities  (i.e.,  the 
ratio  of  the  wave  energy  divided  by  the  ther¬ 
mal  energy  of  the  plasma)  generated  as  a  result 
of  the  various  instabilities  produced  as  the  jet 
penetrates  the  ambient  medium. 

The  principal  plasma  waves  can  be  char¬ 
acterized  as  follows:  the  two  stream  instabil¬ 
ity  waves,  Wu  interact  directly  with  the  ambi¬ 
ent  medium,  and  are  “predated”  (principally) 
by  the  oscillating  two  stream  instability  waves, 
W2,  and  the  ion-acoustic  waves,  W5.  To  reit¬ 
erate,  these  waves  are  generated  by  instabili¬ 
ties  produced  by  the  jet-ambient  medium  in¬ 


teractions.  The  waves  generated  in  the  plasma 
by  the  jet  produce  regions  of  high  electric 
field  strength  and  relatively  low  density,  the 
so-called  “cavitons”  (after  solitons  or  soli¬ 
tary  waves)  which  propagate  like  wave  pack¬ 
ets.  These  cavitons  mix,  collapse,  and  reform, 
depositing  energy  into  the  ambient  medium, 
transferring  momentum  to  it,  and  entraining 
(i.e.,  dragging  along  and  mixing)  the  ambient 
medium  within  the  jet.  The  typical  caviton  size 
when  formed  is  of  order  lO’s  of  Debye  lengths, 
where  a  Debye  length,  Ad  -  7.43x10^ 
cm,  T  is  the  electron  temperature  in  units  of 
eV,  and  Up  is  the  number  density  of  the  ambi¬ 
ent  medium  in  units  of  cnT^ . 

In  order  to  determine  the  energy  deposition 
rate,  the  momentum  transfer  rate,  and  heat¬ 
ing,  we  model  the  plasma  interaction  as  a  sys¬ 
tem  of  very  stiff,  coupled  differential  equations 
(see,  e.g.,  Beall  et  al.  2006),  which  simulate  the 
principal  elements  of  the  plasma  processes  that 
draw  energy  out  of  the  jet.  As  a  test  of  the  fealty 
of  this  method,  we  “benchmark”  (see  Oreskes 
et  al.  1994)  the  wave  population  code  by  using 
the  PIC  code  in  regions  of  the  parameter  space 
where  running  the  PIC  code  simulation  is  prac¬ 
ticable.  We  then  use  the  wave  population  code 
for  regions  of  more  direct  astrophysical  inter¬ 
est.  A  more  detailed  discussion  of  the  compar¬ 
isons  between  the  PIC-code  simulations  and 
the  wave-population  model  can  be  found  in 
Rose,  Guillory,  and  Beall  (2002,  2005).  The 
coupling  of  these  instability  mechanisms  is  ex¬ 
pressed  in  the  model  through  a  set  of  rate  equa¬ 
tions.  These  equations  are  discussed  in  some 
detail  by  Rose  et  al.  (1984,  1987),  and  Beall 
(1990). 

The  solutions  to  the  wave  population  equa¬ 
tions  give  a  normalized  wave  energy.  This 
wave  energy  density  is  then  used  to  estimate 
the  energy  deposition  rate  of  the  jet  into  the 
ambient  medium,  the  propagation  length  of  the 
jet,  the  heating  of  the  plasma,  and  the  momen¬ 
tum  transfer  rate  from  the  jet  to  the  plasma  (see 
Beall  et  al.  2006,  2009  for  a  detailed  discus¬ 
sion). 

Beall  et  al.,  (2006)  illustrates  two  possible 
solutions  for  the  system  of  coupled  differential 
equations  that  model  the  jet-ambient  medium 
interaction:  a  damped  oscillatory  and  an  oscil- 
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latory  solution.  The  Landau  damping  rate  for 
the  two-temperature  thermal  distribution  of  the 
ambient  medium  is  used  for  these  solutions.  As 
noted  in  the  figure  caption,  transitions  toward 
chaotic  solutions  have  been  observed  for  very 
large  growth  rates  for  the  two-stream  instabil¬ 
ity. 

In  order  to  benchmark  the  wave-population 
code,  we  use  that  code  to  calculate  the  propa¬ 
gation  length  of  an  electron-positron  jet  as  de¬ 
scribed  above.  Specifically,  we  model  the  in¬ 
teraction  of  the  relativistic  jet  with  the  ambient 
medium  through  which  it  propagates  by  means 
of  a  set  of  coupled,  differential  equations  which 
describe  the  growth,  saturation,  and  decay  of 
the  three  wave  modes  likely  to  be  produced  by 
the  jet-medium  interaction.  First,  two-stream 
instability  produces  a  plasma  wave,  Wi,  called 
the  resonant  wave,  which  grows  initially  at 
a  rate  Fi  <  ( V3/27)(ni/2n^)'^^Wp,  where 
7  is  the  Lorentz  factor  of  the  beam,  n/,  and 
Hp  are  the  beam  and  cloud  number  densi¬ 
ties,  respectively,  and  aip  is  the  plasma  fre¬ 
quency,  as  described  more  fully  in  Rose  et 
al.  (1984).  dEpiasmaldx  -  -(llnbVb)(daeildt), 
can  be  obtained  by  determining  the  change 
in  7  of  a  factor  of  2  with  the  integration 
j dy  -  -  J[d(aei)/dt]/(vbnhm' c^)  as  shown 
in  Rose  et  al.,  1978  and  Beall  1990,  where 
m  is  the  mass  of  the  beam  particle.  Thus, 
Lp  -  {{1 12)ycnbmc^) I {daei  I dt)  cm  is  the  char¬ 
acteristic  propagation  length  for  collisionless 
losses  for  an  electron  or  electron-positron  jet, 
where  daeijdt  is  the  normalized  energy  depo¬ 
sition  rate  (in  units  of  thermal  energy)  from 
the  plasma  waves  into  the  ambient  plasma.  In 
many  astrophysical  cases,  this  is  the  dominant 
energy  loss  mechanism. 

The  average  energy  deposition  rate,  < 
d(ae\)ldt  >,  of  the  jet  energy  into  the  ambi¬ 
ent  medium  via  plasma  processes  can  be  cal¬ 
culated  as  <  d{aei)ldt  >=  UpkT  <  IT  > 
(ri/ajp)a)p  ergs  cmT^ s^^ ,  where  rip  is  num¬ 
ber  density  in  units  of  cm-3  of  the  ambient 
medium,  k  is  Boltzmann’s  constant,  T  is  the 
plasma  temperature,  <  IT  >  is  the  average  (or 
equilibrium)  normalized  wave  energy  density 
obtained  from  the  wave  population  code,  F i  is 
the  initial  growth  rate  of  the  two-stream  insta¬ 
bility,  and  ojp  is  the  plasma  frequency. 


It  is  of  interest  to  compare  the  results  of 
a  Particle-In-Cell  (PIC)  code  simulation  of  an 
electron-positron  jet  propagating  through  an 
ambient  medium  of  an  electron-proton  plasma 
with  the  solutions  obtained  by  the  wave  pop¬ 
ulation  model  code.  A  small  magnetic  field  is 
applied  along  the  jet’s  longitudinal  axis  to  sup¬ 
press  a  filamentation  instability,  but  this  does 
not  affect  the  propagation  length,  Lp,  which  is 
our  principal  concern  here.  Lp  is  the  distance 
over  which  the  plasma  instabilities  reduce  the 
beam  gamma  by  a  factor  of  two.  At  the  same 
time,  the  ambient  medium  is  heated  and  en¬ 
trained  into  the  jet.  We  believe  that  this  config¬ 
uration  is  a  reasonable  end  point  for  the  initial 
interaction  of  the  relativistic  jet  with  the  inter¬ 
stellar  medium,  given  the  pressure  exerted  on 
the  ambient  medium  with  an  oblique  or  trans¬ 
verse  magnetic  field.  These  simulations  show 
that  a  relativistic,  low-density  jet  can  interpen¬ 
etrate  an  ambient  gas  or  plasma. 

Figure  1  shows  the  wave  levels  as  cal¬ 
culated  dynamically  from  our  solution  to  the 
wave  population  code  as  compared  to  the  PIC- 
code  simulation.  The  plasma  density  for  this 
simulation  is  rip  -  IcmT^,  the  plasma  temper¬ 
ature  is  Tc  -  10‘^K,  the  ratio  of  the  beam  den¬ 
sity  to  the  plasma  density,  R  -  10“^,  and  the 
plasma  has  a  hot  electron  tail  (produced  by  the 
jet)  with  a  temperature  of  10^  K.  The  vertical 
axis  is  the  normalized  wave  energy  density,  and 
the  horizontal  axis  is  time,  expressed  in  units 
of  plasma  periods,  where  the  plasma  period, 
cjp  -  5.64x10'^  ynjj. 

Initially,  and  for  a  significant  fraction  of  its 
propagation  length,  the  principal  energy  loss 
mechanisms  for  such  a  jet  interacting  with  the 
ambient  medium  is  via  plasma  processes  (Rose 
et  al.  1984,  Beall  1990). 

3.  PIC  code  simulations  showing 
momentum  transfer  to  the  ambient 
medium  from  beam-generated 
plasma  waves 

As  part  of  our  research  into  the  micro-physics 
of  the  interaction  of  jets  with  an  ambient 
medium,  we  have  begun  to  investigate  in  some 
detail  the  transfer  of  momentum  from  the  jet. 
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Jet25  —  y  =2  election-position  beams,  n^/^  =  20, 


Fig.  1.  This  figure  shows  a  recent  PIC  simulation  of  a  jet-ambient-medium  with  the  propagation  length  for 
the  wave  population  code  marked,  and  for  7  =  2,  =  20,  and  ~  4Bc.  As  can  be  seen,  the  propagation 

length,  Lp,  calculated  in  Beall  et  al.,  2005,  is  in  agreement  with  the  station  along  the  axis  of  the  jet  where 
the  plasma  waves  have  developed  significantly. 


Understanding  how  such  a  transfer  is  accom¬ 
plished  is  essential  to  understanding  the  man¬ 
ner  in  which  the  ambient  medium  (for  exam¬ 
ple,  from  interstellar  clouds)  is  accelerated  and 
eventually  entrained  into  the  large-scale  astro- 
physical  jet. 

We  intend  a  further  analysis  of  the  PIC 
code  simulations  presented  here  in  order  to 
study  several  issues  related  to  caviton  dynam¬ 
ics,  including  momentum  transfer  from  the 
cavitons  to  the  background  thermal  plasma, 
and  energetic  electron  and  proton  acceleration. 

As  a  further  continuation  of  this  research, 
we  plan  to  develop  a  multi-scale  code  which 
uses  the  energy  deposition  rates  and  momen¬ 
tum  transfer  rates  from  the  PIC  (Particle-In- 
Cell)  and  wave-population  models  as  source 
terms  for  the  highly  parallelized  hydrodynamic 
code  currently  running  on  the  NRL  Cray  XD- 
1. 


We  use  a  PIC  code  simulation  to  model 
caviton  formation  and  momentum  transfer. 
This  simulation  is  a  one-dimensional  (1-D) 
simulation  with  a  length,  L  -  6x10®  cm,  us¬ 
ing  periodic  boundaries.  In  the  simulation,  L  is 
greater  than  the  maximum  two-stream  unstable 
wavelength,  the  ratio  of  plasma  to  beam  densi¬ 
ties  =  ripeltihe  -  20,  Hp  -  Icm  -  3,  the  beam 
velocity,  Vb  =  0.866c,  and  the  grid  resolution 
0.2xAd,  the  Debye-length,  in  order  to  resolve 
the  caviton  structures.  It  is  an  understatement 
of  considerable  magnitude  to  say  that  this  cal¬ 
culation  is  computationally  intensive. 

This  result  shows  clearly  the  acceleration 
of  the  ambient  medium  via  plasma  processes 
in  the  interaction  of  an  astrophysical  jet  with 
the  ambient  medium.  It  is  the  first  such  result 
that  explicitly  demonstrates  momentum  trans¬ 
fer  from  a  consideration  of  the  micro-physics 
of  the  interaction.  We  intend  to  use  these  re- 
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Sample  c^itoii  formation  illustrates  net  positive 
axial  momentum  of  density  structure: 


I 

a 


Sample  caviton : 


Ion  phase  space  contours 
show  that  caviton  has  a  netforward 
directed  momentum : 


Fig.  2.  PIC  code  simulation  which  for  the  first  time  specifically  identifies  one  of  the  cavitons  formed  by 
a  jet-ambient  medium  interaction.  The  left  panel  shows  the  electric  field  intensity  (in  kV/cm)  along  the 
x-axis.  The  right  panel  plots  the  velocity  distribution  of  ions  in  the  ambient  medium  vertically  as  a  function 
of  velocity,  and  as  a  function  of  position  along  the  jet  axis  horizontally.  Note  the  white  region  at  roughly 
5x10^  cm  indicates  a  deficiency  in  the  density  of  the  ambient  medium  about  a  zero  velocity.  This  indicates 
that  the  caviton  at  the  same  place  in  the  jet  has  given  the  ambient  medium  a  positive  velocity  at  that  location 
in  the  jet. 


suits  to  aid  us  in  benchmarking  an  analytical 
calculation  of  the  momentum  transfer  so  that 
we  can  incorporate  the  plasma  processes  as 
a  source  term  for  the  hydrodynamic  code  we 
have  ported  to  the  NRL  Cray  XD-1. 

4.  PIC  code  modeling  of  the 
interaction  of  a  plasmoid  with  a 
fixed  magnetic  field 

It  is  possible  (using  PIC  code  simulations 
on  available  computing  platforms)  to  investi¬ 
gate  how  jets  might  interact  with  magnetized 
clouds.  Figure  3  shows  a  frame  from  a  “movie” 
made  from  a  PIC  code  simulation  for  a  related 
problem  of  an  electron-proton  plasmoid  inter¬ 
acting  with  a  fixed  magnetic  field.  P.  The  par¬ 
ticular  frame  shown  in  the  figure  demonstrated 


the  significant  penetration  of  the  magnetic  field 
my  the  plasma,  and  shows  that  the  interaction 
actually  accelerated  some  of  the  plasmoid’s 
particles. 

We  are  considering  the  possibility  of  addi¬ 
tional  simulations  of  this  sort  as  a  means  of 
modeling  cross-field  injection  of  jets  into  the 
magnetic  fields  of  interstellar  clouds.  But  even 
at  first  gloss,  it  is  apparent  that  the  plasmoid 
can  interpenetrate  the  even  a  strong  magnetic 
field,  and  that  the  interaction,  in  the  center  of 
momentum  frame  of  the  system  of  a  jet  and  an 
interstellar  cloud,  can  accelerate  the  cloud.  We 
believe  that  the  outcome  of  such  interactions 
will  be  a  magnetic  field  in  the  interstellar  cloud 
co-aligned  with  the  jet. 

It  is  not  too  much  to  say  that  simulations  of 
such  interactions  using  PIC  codes  is  computa¬ 
tionally  intensive. 
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Fig.  3.  Illustration  of  the  capabilities  of  a  PIC  code  simulation  to  model  the  interaction  of  an  electron- 
proton  plasmoid  with  a  fixed  magnetic  field.  The  simulation  shows  the  penetration  of  the  magnetic  field  by 
the  plasmoid.  Magnetic  fields  “anchored”  in  an  interstellar  cloud  are  thus  likely  to  entrain  and  accelerate 
the  cloud  as  they  penetrate  the  magnetic  field  and  force  it  into  an  orientation  aligned  with  the  jet. 


5.  Hydrodynamic  modeling  of 
supersonic  jet  in  a  Seyfert  galaxy 

As  part  of  our  project  to  incorporate  the  plasma 
effects  into  a  multi-scale  code,  we  have  be¬ 
gun  modeling  the  interaction  of  a  supersonic 
jet  (for  example,  a  Seyfert  galaxy  jet),  using  a 
highly  parallelized  version  of  the  VH-1  code 
(see  Figure  4). 

The  code  has  been  ported  to  the  NRL  Cray 
XD-1  as  part  of  a  collaboration  with  Curtis 
Saxton  and  Kinwah  Wu  at  the  Mullard  Space 
Sciences  Laboratory,  in  order  to  model  the  in¬ 
teraction  of  a  Seyfert  jet  with  the  interstellar 
medium.  We  have  begun  initial  runs  at  various 
resolutions  up  to  a  3-dimensional  simulation 
of  1024^  cells.  This  simulation  will  involve  a 
calculation  of  more  than  IOjcIO®  floating  point 
numbers  for  each  time  ste. 


6.  Concluding  remarks 

Plasma  effects  can  have  observational  conse¬ 
quences.  Beall  (1990)  has  noted  that  plasma 
processes  can  slow  the  jets  rapidly,  and  Beall 
and  Bednarek  (1999)  have  shown  that  these  ef¬ 
fects  can  truncate  the  low-energy  portion  of  the 
y-rays  spectrum  (see  Figure  3).  In  the  inter¬ 
ests  of  brevity,  we  do  not  go  into  the  (reason¬ 
able)  assumptions  for  the  calculation.  Please 
see  Beall  and  Bednarek  (1999)  for  a  detailed 
discussion.  A  similar  effect  will  occur  for  neu¬ 
trinos  and  could  also  reduce  the  expected  neu¬ 
trino  flux  from  AGN.  The  presence  of  plasma 
processes  in  jets  can  also  greatly  enhance  line 
species  by  generating  high-energy  tails  on  the 
Maxwell-Boltzmann  distribution  of  the  ambi- 
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Paiallelized  VH-1  Laige -Scale  Jet  Simulation 
shown  witli  a  volumehic  rendering 


j  et  density 


jet  presstjre 


Fig.  4.  Volumetric  rendering  of  the  density  and  pressure  of  the  initial  stages  of  jet  in  a  Seyfert  galaxy,  using 
the  highly-parallized  version  of  the  VH-1  code  (in  collaboration  with  Curtis  Saxton  and  Kinwah  Wu).  The 
images  show  only  the  outermost  shell  of  the  ambient  medium  as  it  expands  into  the  surrounding  ISM.  The 
actual  jet  structure  is  buried  in  the  core  of  this  volume. 


ent  medium,  thus  abrogating  the  assumption 
of  thermal  equilibrium.  Beall,  Guillory,  and 
Rose  (1999)  have  also  calculated  the  greatly 
enhanced  line  emission  for  certain  species  in 
the  presence  of  a  relativistic  jet  due  to  plasma 
processes. 

The  hypothesis  of  jets  from  AGN  interact¬ 
ing  with  the  intracluster  medium  via  collision¬ 
less  (plasma)  processes  requires  that  the  jets 
overcome  collisional  and  collisionless  losses 
and  propagate  to  significant  distances  into  the 
intracluster  medium.  This  in  turn  allows  us  to 
constrain  the  jet  parameters  as  the  jet  emerges 
from  the  elliptical  AGN.  In  general,  the  jets 
must  have  values  of  y,  the  ratio  of  the  total 
particle  energy  over  the  particle  rest  mass,  that 
are  at  least  rather  relativistic  over  a  significant 
fraction  of  their  propagation  length. 

An  analysis  of  the  energy  loss  due  to 
plasma  processes,  taken  from  the  preceding 


equations  and  references,  and  the  computer 
simulations  that  determine  (daeildt),  the  av¬ 
erage  wave  energy  deposition  into  the  ambi¬ 
ent  medium  per  unit  time,  yields  some  use¬ 
ful  bounds  for  possible  energy  deposition  rates 
due  to  plasma  processes.  We  can  further  con¬ 
strain  the  jet  parameters  by  expressing  the  ki¬ 
netic  luminosity  of  the  jet  as  Ph  =  dE/dt  - 
ymc^rihVhnr^,  where  y  is  the  ratio  of  total  en¬ 
ergy  to  rest  mass  energy,  mc^  is  the  rest  mass 
energy  of  the  beam  particles,  v*  is  the  beam  ve¬ 
locity,  and  is  the  beam  radius. 

If  the  beam  is  significantly  heated  by  the 
jet-cloud  interaction,  the  beam  will  expand 
transversely  as  it  propagates,  and  will  there¬ 
fore  have  a  finite  opening  angle.  These  “warm 
beams”  result  in  different  growth  rates  for 
the  plasma  instabilities,  and  therefore  pro¬ 
duce  somewhat  different  propagation  lengths 
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(see,  e.g.,  Kaplan  and  Tsytovitch,  1973,  Rose, 
Guillory,  and  Beall,  2002,  and  Beall  et  al. 
2006).  A  “cold  beam”  is  assumed  to  have  lit¬ 
tle  spread  in  momentum.  The  likely  scenario 
is  that  the  beam  starts  out  as  a  cold  beam  and 
evolves  into  a  warm  beam  as  it  propagates 
through  the  ambient  medium.  This  scenario  is 
clearly  illustrated  by  the  Particle-In-Cell  (PIC) 
simulations  we  have  used  to  benchmark  the 
wave  population  codes  appropriate  for  the  as- 
trophysical  parameter  range  (see,  e.g.,  Beall, 
Guillory,  and  Rose,  1999,  and  Rose,  Guillory, 
and  Beall,  2005). 

Assuming  that  the  ambient  medium  is  also 
significantly  heated  by  the  jet  (at  some  late 
times  in  the  history  of  the  interaction),  the  am¬ 
bient  medium  will  develop  in  effect  a  two  tem¬ 
perature  distribution  (see,  e.g.,  Beall,  Guillory, 
and  Rose,  1999).  The  effect  of  the  jet-cloud  in¬ 
teraction  is  to  produce  a  high-energy  tail  on  the 
thermal  distribution  of  electron  energies  in  the 
cloud.  This  high-energy  tail  critically  alters  the 
Landau  damping  rate  of  waves  in  the  plasma. 
An  analytical  calculation  of  the  boost  in  en¬ 
ergy  of  the  electrons  in  the  ambient  medium 
to  produce  such  a  high  energy  tail,  with  Eket  ~ 
30  -  lOOkT,  is  confirmed  by  PIC -code  simula¬ 
tions.  Aside  from  altering  the  Landau  damp¬ 
ing  rate,  such  a  high-energy  tail  can  greatly 
enhance  line  radiation  over  that  expected  for 
a  thermal  equilibrium  calculation  (see  Beall  et 
al.  2006,  and  Beall,  Guillory,  and  Rose  1999 
for  a  detailed  discussion). 

The  dominant  energy  loss  mechanism  for 
the  jets  propagating  outward  from  the  AGN 
in  the  cluster  core  will  be  due  to  plasma  pro¬ 
cesses  for  some  significant  portion  of  the  jet’s 
propagation  length,  and  before  the  jet  transi¬ 
tions  to  a  hydronamic  regime.  Electron-hadron 
jets  can  apparently  propagate  to  distances  of 
10s  to  100s  of  Kpc,  despite  plasma  processes 
being  dominant  over  other  energy  loss  mech¬ 
anisms  in  most  cases.  In  an  electron-proton 
(or  electron-hadron)  jet,  the  electrons  lose  en¬ 
ergy  to  plasma  processes  more  rapidly  than  do 
the  protons.  The  jet  protons  therefore  drag  the 
electrons.  This  produces  a  current  along  the  jet 
in  the  jet’s  rest  frame.  A  magnetic  field  so  pro¬ 
duced  will  stabilize  the  jet. 

The  presence  of  hadrons  in  the  jet  will  pro¬ 
duce  nuclear  y-rays  and  neutrinos  as  it  inter¬ 


acts  with  the  ambient  medium  (see  Beall  and 
Bednarek,  1999  for  a  discussion).  The  plasma 
instabilities  modify  the  emitted  y-ray  spectrum 
significantly.  If  jets  are  hadronic  (a  scenario 
that  would  help  with  both  the  energy  transport 
problem  and  their  propagation  length),  then 
they  probably  also  have  a  significant 
component  that  will  ‘fill  in”  to  account  for 
some  of  the  observed  radiation. 

The  detection  of  neutrinos  from  jet  sources 
would  directly  suggest  an  hadronic  component 
at  relativistic  (i.e.  early)  stages  of  jet  formation. 
Such  an  hypothesis  is  consistent  with  Eichler’s 
(1979)  suggestion  of  using  neutrinos  as  a  probe 
of  AGN.  The  recent  results  from  the  Auger 
Collaboration  (Letessier-Salvon  2009)  suggest 
an  extremely  energetic  hadronic  component  as 
a  consequence  of  these  jets. 

7.  DISCUSSION 

SERGIO  COLAFRANCESCO:  What  is  the 
efficiency  of  the  formation  of  the  cavitons  on 
the  jet? 

JIM  BEALL:  The  cavitons  are  ubiquitous 
during  the  jet  interaction.  They  are  continually 
formed  (being  generated  from  different  wave 
length  plasma  waves),  phase  mix  (because  the 
plasma  is  a  dispersive  medium),  and  then  col¬ 
lapse.  The  process  of  collapse  is  likely  to  be 
much  like  waves  breaking  on  a  shore,  so  they 
are  likely  to  transfer  significant  momentum  to 
the  ambient  medium.  But  we  need  to  model 
this  quantitatively  to  be  certain  of  how  much 
momentum  is  transferred. 

DMITRI  BISIKALO:  What  happens  with  the 
magnetic  field  in  the  cloud  during  its  interac¬ 
tion  with  the  jet? 

JIM  BEALL:  It’s  clear  from  the  movie  of 
our  simulation  (which  was  done  for  another 
project,  by  the  way)  that  the  interaction  is  in¬ 
credibly  violent.  It  is  plausible  that  the  jet 
transfers  momentum  to  the  magnetic  field, 
and  tends  to  force  the  transverse  magnetic 
fields  to  co-align  themselves  with  the  longitu¬ 
dinal  axis  of  the  jet.  This  general  picture  seems 
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right  to  me,  but  the  details  are  very  difficult  to 
model. 

SERGIO  COLAFRANCESCO:  What  is  the 
effect  of  magnetic  reconnection  on  the  forma¬ 
tion  of  the  cavitons  in  the  jet? 

JIM  BEALL:  The  cavitons  are  of  order  cen¬ 
timeters  in  size,  so  they  will  form  without  in¬ 
terference  from  the  magnetic  fields.  However, 
other  instabilities  are  affected  by  the  magnetic 
field.  For  example,  a  longitudinal  magnetic 
field,  like  the  ones  inferred  from  polarization 
measurements,  tends  to  suppress  a  filamenta- 
tion  instability  that  causes  the  jet  to  thread  into 
many  more,  smaller-diameter  jets  which  will 
then  can  recombine.  This  process  does  not  af¬ 
fect  estimations  of  the  energy  deposition  rate 
or  the  propagation  length  of  the  jet. 
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